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ABSTRACT 

We study the temporal evolution of umbral dots (UDs) using measurements from the CRISP imaging 
spectropolarimeter at the Swedish 1-m Solar Telescope. Scans of the magnetically sensitive 630 nm 
iron lines were performed under stable atmospheric conditions for 71 min with a cadence of 63 s. 
These observations allow us to investigate the magnetic field and velocity in and around UDs at 
a resolution approaching 0^13. From the analysis of 339 UDs, we draw the following conclusions: 
(1) UDs show clear hints of upfiows, as predicted by magnetohydrodynamic (MHD) simulations. By 
contrast, we could not find systematic downfiow signals. Only in very deep layers we detect localized 
downfiows around UDs, but they do not persist in time. (2) We confirm that UDs exhibit weaker and 
more inclined fields than their surroundings, as reported previously. However, UDs that have strong 
fields above 2000 G or are in the decay phase show enhanced and more vertical fields. (3) There are 
enhanced fields at the migration front of UDs detached from penumbral grains, as if their motion were 
impeded by the ambient field. (4) Long-lived UDs travel longer distances with slower proper motions. 
Our results appear to confirm some aspects of recent numerical simulations of magnetoconvection in 
the umbra (e.g., the existence of upfiows in UDs), but not others (e.g., the systematic weakening of 
the magnetic field at the position of UDs.) 
Subject headings: Sun: magnetic fields - sunspots - convection 



1. INTRODUCTION 

Umbral dots (UDs) are transient brightenings observed 
in sunspot umbrae and pores, with typ ical sizes of 300 km 
and lifetimes of 10 min (e.g., (Sobotka et al. 1997a|b ). 
They cover only 3-10% of the umbral area, but con- 
tribute 10-20% of its brightness. For this reason, UDs 
are believed to play a vital role in the energy balance 



of sunspots (Deinzer 1965 Sobotka et al. 1993 
etlI][2QTQ] ). 
UDs exhibi 



Moradi 



ibit systematic proper motions in mature 
sunspots: those appearing in the central umbral region 
are static, while UDs in peripheral reg ions move inward 
with an average velocity of 1.0 kms~^ ( Kitai et aL||2007 
Riethmiiller et ah] 2008b). Some peripheral UDs are 
the continuation of penumbral grains — bright elongated 
structures at the head of penumbral filaments that move 
toward the center of the sunspot with speeds of about 
0.4kms~^ dSobotka et al. ' ' 



1999a Rimmele & Marino 



2006). When the migration front detaches into a circular 



bright point, the tip of the penumbral grain becomes an 
UD. 

It is believed that the mechanism behind UDs is con- 
vection interacting with the strong vertical field of the 
umbra, and many observational results support this idea 



([Riethmuller et al.l|2008a[ [Bharti et al 
abe et al. ' 2010p ! In the formation pha^ 



|2010[ |Watan- 

lase of sunspots, 

UDs are akin to granules but their apparent motion is 
more stochastic because of the surrounding magnetic 
field ( Sobotka et al.|1999b ). In developed sunspots, UDs 
are small and quiescent due to the stronger suppression 
of convection. 

UD research is entering a new phase in which computer 
simulations gu ide observational efforts. The innovative 
2006) predicted UDs 

with central 



localized downfiow 



k lanes and small 
patches at their ends. A clear detection of those fea- 
tures would immediately validate the numerical models, 
so they have been the target of recent observations. The 
dark lanes are the result of enhanced density in the up- 
per central part of UDs, caused by the piling up of hot 
gas that rises from deeper down. Once the gas reaches 
the surface, it cools by radiative losses and descends in 
narrow downfiow channels at the end of the dark lanes. 



Bharti et al. (2007b) observed a dark lane in a big UD. 
However, the very large size of this UD (> 1000 km) sug- 
gests that it could actually have been a cluster of sev- 
eral UDs. The first detec tion of downfiows sur rounding 
an UD was presented by Bharti et al. (2007a). Subse- 
quently, Ortiz et al. ( 2010p reported solid evidence of 
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dark lanes and localized downflows based on spectropo- 
larimetric observations taken at the Swedish 1-m Solar 
Telescope. The s izes of the dark la nes and downflowing 
patches found by Ortiz et al. (2010) are near the diffrac- 
tion limit of the telescope, with the substructures keeping 
their identity for periods of only a few minutes. These 
authors also reported enhanced net circular polarization 
at the site of the downflows. 

The evolution of UDs and their magnetic fields is dif- 
ficult to study — and hence poorly known — because one 
needs full vector spectropolarimetric measurements at 
very high temporal and spatial resolution. To the best 
of our knowledge, the magnetic properties of UDs have 
never been investigated at the required cadence a n d spa- 
tial r esolution (but see Sobotka & Jurcak'"2009'). Ortiz 
et al.| (2010) performed a preliminary analysis of the tem- 
poral evolution of six UDs, and this work should be con- 
sidered a substantial extension of their study. Both the 
cadence and the polarimetric sensitivity of o ur measure- 
ments are improved with respect to those of [Ortiz et al.| 
(^2010), as is the total duration of the observations, 71 
minutes, during which the seeing conditions were excel- 
lent and stable. We use this unique data set to investigate 
the evolution of the magnetic and velocity fields in and 
around UDs. 

The paper is organized as follows. The observations 
are described in Section [2j followed by an account of the 
methods used for the detection of UDs and derivation 
of the velocity and magnetic information (Section [3|. In 
Section [4] we quantify how convection is modified m the 
umbra. In Section |5] we describe the evolution of some 
typical UDs, and present the results of our statistical 
analysis. Finally, based on these results, we discuss the 
physical properties of UDs in Section [6j 

2. OBSERVATIONS AND DATA REDUCTION 

The observations were obtained with the CRisp Imag- 
ing Spectro-Polarim eter (CRISP) at the S wedish 1-m So- 
lar Telescope (SST; [Scharmer et aL]|2003^ on La Palma, 
Spain. CRISP is based on a dual Fabry-P^rot interferom- 
eter similar to that described by Sch armer, (200 6). The 
incoming light is modulated by two liquid crystiQ variable 
retarders cycling through four states and then analyzed 
by a polarizing beam splitter in front of two narrow-band 
cameras. The narrow-band cameras record orthogonal 
polarization states to minimize seeing-induced crosstalk. 
CRISP has a third camera for wide band imaging. All 
the cameras operate at 35 frames per second and take 
exposures of 17.6 ms. The synchronism between them is 
ensured by an external optical chopper. 

CRISP was used to measure the four Stokes profiles 
of the magnetically sensitive Fe I 6301.5 and 6302.5 A 
lines, each sampled at 15 wavelength positions in steps of 
48 mA, from —350 to +322 mA. Line positions 0-14 sam- 
ple the 6301.5 A line, while positions 15-29 correspond 
to the 6302.5 A line. In addition, a continuum wave- 
length point (6303.2 A, line position 30) was measured. 
We recorded 9 frames per modulation state, resulting in 
36 exposures per wavelength position. The total time for 
a complete wavelength scan of the two Fe I lines plus the 
continuum point was 32 s. Another 30 s were needed to 
scan the Ca II line at 8542 A (not used in this paper). 
Thus, the temporal cadence of the Fe I scans is 63 s. The 
average CRISP transmission profile has a Gaussian core 



(FWHM of 64 mA at 6300A) and wide Lorentzian wings. 
This transmission profile reduces the line core depth of 
the Fe I lines by about 20%. 

The CRISP etalons are mounted in tandem on a tele- 
centric beam. The separation of the cavities of the high 
resolution etalon sets the wavelength of the transmis- 
sion profile. This is not strictly the same over the whole 
field of view (FOV) because the surface cannot be in- 
finitely flat, producing random wavelength shifts (cav- 
ity errors) across the FOV. Intensity fluctuations, intro- 
duced by cavity errors in the presence of a spectral line, 
were removed fr om the fiat-fie l d ima ges in a similar way 
as described by Schnerr et al. (2011|). Their flat-fielding 
scheme computes a cavity-error-free averaged quiet-sun 
profile, which is removed from the flats on a pixel- by- 
pixel basis. This average is obtained by summing many 
hundreds of exposures acquired by moving the telescope 
in circles around the disk center. The difference between 



the Schnerr et al. (2011) scheme and ours is that we ap- 
plied a polarization-free fiat for all the four Stokes states 
taken at one wavelength position. 

We also corrected the data for spectral intensity gradi- 
ents introduced by the CRISP prefilter. The prefilter cor- 
rection can be decomposed in two contributions, an av- 
erage prefilter shape and a term that accounts for pixel- 
to-pixel deviations: 



P(A,x,^) = Pg(A)+dP(x,^). 



(1) 



The second term of the equation is included in the fiat- 
field correction. To obtain the global prefilter shape 
(Pg), our estimate of the averaged quiet-sun profile is 
compared to the FTS atlas convolved with the CRISP 
transmission profile. The prefilter is assumed to have a 
Lorentzian shape which is multiplied with a polynomial 
term to account for asymmetries: 



^g(A) = 



1 + [2(A- Ao)/H''^- 



(1+PoA+PiA^). (2) 



Here, Aq and w are the central wavelength and FWHM 
of the prefilter, respectively, A^cav is the number of cav- 
ities of the filter (A/'cav = 2 in this case), and po and pi 
are the coeffici ents of the polyno mial. A least-squares- 
fitting scheme ( [Markwardt 2009) was used to compute 
the prefilter parameters, minimizing the differences be- 
tween the observed and modeled curves. All Stokes pa- 
rameters were then corrected by dividing each spectrum 
with the fitted prefilter curve. 

The theoretical diffraction limit of the telescope around 
6300 A is 0^13, and the image pixel size is 0^.^059. To en- 
sure very high spatial resolution, we used the adaptive 
optics system of the SST and the Multi-Object M ulti- 
Frame Blind Deco nvolution technique (MOMFBD; |van| 
Noort et al]|2005| . The MOMFBD algorithm considers 
all frames taken in one scan (31 wavelength points x 4 
modulation states x 9 repetitions x 3 cameras) to re- 
move image distortions from the individual filtergrams. 

Under enhanced differential-seeing conditions, residual 
rubber sheet distortions are present between the narrow- 
band images of the same scan even after MOMFBD 
restoration. This effect appears when the size of the 
patches used to divide the images for the MOMFBD 
processing (here 128 x 128 pixels) is larger than the spa- 
tial scale of the atmospheric distortion (or iso-planatic 
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Fig. 1. — Filtergrams from the best scan of the data set, taken at 08:30 UT. The full FOV is shown. From left to right: Stokes / in 
the blue wing of 6301.5 A (line position 0), Stokes / at line center (line position 7), and Stokes V in the blue wing (line position 5). The 
direction to disk center (DC) is displayed with an arrow. The white rectangle indicates the FOV of Figure |2] The bottom row shows the 
Stokes profiles emerging from the UD marked with triangles in the upper panels. The +-symbols indicate the measured signals. 

patch). To attenuate their influence in our measure- 
ments, we employed an extra step in the processing fol- 
lowing an idea from V. Henriques (private communi- 
cation). The wide-band images were used twice in a 
MOMFBD restoration in the following manner: 

1. All the frames were combined to produce the ref- 
erence anchor image. 

2. The frames were separated in sets associated with 
each wavelength position and modulation state, re- 
sulting in one restored wide-band image per wave- 
length and modulation state. 

The second set of wide-band images was not used 
for wavefront sensing, i.e., they did not contribute to 
the determination of the Point Spread Function in the 
MOMFBD calculations. The state-dependent restored 
wide-band images were compared with the anchor wide- 
band image to remove the residual rubber sheet deforma- 
tions in the individual flltergrams. This correction was 
applied prior to the demodulation of the data to achieve 
almost perfect co-alignment between the 4 modulation 
states from which the Stokes parameters are derived at 
any wavelength position. The method also achieves al- 
most perfect co-alignment between the different wave- 
length positions to ensure the integrity of the Stokes pro- 
flies over the FOV. 

After restoration, the images were demodulated and 
corrected for instrumental polarization u sing the tele- 
scope model developed by Selbing (2010). For details. 



see van Noort & Rouppe van der VoortN2008fc. In addi- 
tion, we corrected small residual crosstalks from Stokes 
/ to Q, and V by forcing the polarization to be zero 
in the far line wings. All the Stokes proflles were normal- 
ized to the average quiet-sun continuum intensity (/qg) 
computed for each time step. The typical noise levels 
in Stokes Q, U, and V are 1.9 x 10"^, 2.8 x 10"^, and 
1.9 X 10~^ of the continuum intensity, respectively. 

We followed the main sunspot of NOAA active region 
11024 from 08:05 to 09:16 UT on 6 July 2009. The spot 
was located at 25° S and 23° W (heliocentric angle of 36°, 
/i=0.81). During the observations the atmospheric con- 
ditions were excellent and stable. Figure [l] shows selected 
flltergrams from one of the best scans (08:30 UT) and the 
four Stokes proflles measured at the position of an UD. 
The size of the fuh FOV is 58'' x 57''. 

Figure [2] displays a close-up of the umbra and the inner 
penumbra on the disk-center side of the spot. The maps 
of mean linear polarization (LP) and circular polarization 
(CP) degree were calculated as 



LP 



/[Q^(A) + [/^(A)]VV/(A)dA 



CP 



JdA 

/inA)|//(A)dA 
JdX 



(3) 



(4) 



with the integration extending over the Fe I 6301.5 A 
line. The LP signals are stronger toward the limb due to 
projection effects. Most of the flne-scale structures seen 
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Fig. 2. — Close-up of the region indicated with a white box in Figure^ Clockwise, starting from upper left: intensity in the blue wing 
of Fe I 6301.5 A (line position 0), continuum intensity at 6303.2 A (line position 30), CP map, and LP map. The direction to disk center 
(DC) is displayed with an arrow. A cut of a point-like feature (corresponding to an UD) is shown in the CP map. The length of the arrow 
in the plot is 0(^13 (equal to the diffraction limit). 

3. The smallest features visible in the maps (UDs and 
dark cores of penumbral filaments) are only slightly 
larger than the diffraction limit of the SST. They 
show higher contrasts in circular polarization. 



in the LP and CP maps coincide with UDs. Figure [2] 
clearly demonstrates that 

1. The first and last intensity images of the scan are 
perfectly aligned and do not show distortions in 
spite of the fact that they were taken about 30 s 
apart (top row). This confirms the good perfor- 
mance of the MOMFBD processing. 

2. Almost all UDs in the umbra keep their brightness, 
size, and position during the line scan (top row), 
that is, the scan time is shorter than the dynamic 
timescale of UDs. 



We believe this is one of the best UD data sets ever ob- 
tained, because of its superb spatial resolution (0^.^13), 
the stable seeing conditions, and the availability of two- 
dimensional spectropolarimetric measurements for more 
than 70 minutes with good temporal cadence (63 s). The 
non-simultaneous acquisition of spectral information — 
the main disadvantage of Fabry-Perot systems — can be 
well neglected because UDs evolve on time scales longer 
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Fig. 3. — Physical parameters derived from the observed Stokes profiles. Shown in the figure is the scan with best seeing conditions, taken 
at 08:30 UT. The top panels, from left to right, display the continuum intensity map with the position of appearance of all the detected UDs 
(yellow circles are central UDs, red squares are peripheral UDs, blue triangles are grain-origin UDs), and two bisector maps sampling high 
('i^high) and deep ('Udeep) photospheric layers. Negative velocities (blue) represent upfiows along the line of sight, while positive velocities 
(red) mean downfiows. The bottom panels, from left to right, show the field strength, the field inclination, and the field azimuth in the 
local reference frame (LRF). The black contours indicate the approximate umbral boundary defined by the spatially-smoothed continuum 
intensity. The arrows mark the direction to disk center (DC). An mpeg version of this figure is available in the on-line journal. 



than the scan time. 

3. DATA ANALYSIS 

In this section we describe the methods we have used 
to detect and track UDs, as well as the line bisector cal- 
culations and Stokes inversions performed to derive their 
velocities and magnetic fields. 

3.1. Detection and Categorization 

For the statistical study of UDs it is convenient to use 
automatic detectio n algorithms (Sobotka et al.^^l997a' 



Bharti et al. 2010). However, we have implemented a 



manual procedure here because, even under the very sta- 
ble seeing conditions of our observations, the image qual- 
ity still shows an unavoidable amount of residual fiuctu- 
ation which might compromise the performance of auto- 
matic methods. 

The procedure works as follows. First we inspect the 
continuum movie to identify the frame of appearance of 
each UD. The position of the UD is then tracked by click- 
ing on the screen until it disappears. After going through 
the temporal sequence, we run the movie backward in 
time to check the consistency. When an UD splits in two 



6 



H. Watanabe et al. 



TABLE 1 
Average properties of UDs 



15 



Parameter 


Central 


Peripheral 


Grain-origin 



^ 


Number 


98 


112 


129 


o 
o 


Lifetime-*^ [min] 


19 


18 


17 


o 


Proper motion [kms"-*^] 


0.19 


0.31 


0.49 




Brightness ratio 


1.46 


1.50 


1.89 




Diameter [km] 


401 


391 


482 





^Excluding UDs which do not appear or disappear within the 
observation period. 

fragments, the bigger coraponent is the one that contin- 
ues to be tracked and the smaher component is selected 
separately as a new entity. When two UDs merge, the 
smaller is assumed to die. In some cases, UDs show re- 
currence at the same position. If the recurrent UDs ap- 
pear within an interval of 3 frames (~3min), we consider 
them a single entity. 

UDs evolving from penumbral grains are also studied. 
In this case, the detection starts from the frame in which 
the tip of the penumbral grain detaches from the fila- 
mentary structure. Sometimes these UDs are connected 
to the penumbral grains through a faint tail, but they 
become more isolated and roundish as they move toward 
the umbra. 

This method has allowed us to obtain the trajecto- 
ries of 339 UDs. According to their place of birth, we 
categorize them into three groups. UDs located in the 
central part of the umbra are called "central UDs" (98 
samples out of 339 UDs), UDs located in the periph- 
eral area are called "peripheral UDs" (112 samples), and 
UDs detached from penumbral grains are called "grain- 
origin UDs" (129 samples). The peripheral area is an 
annular region adjacent to the umbra-penumbra 
boundary. It is known that central UDs are static while 
peripheral and grain-origin UDs show a systeni atic mo- 



tion toward the center of the umbra ( Ewell 1992j Sobotka 
et al.|ll997bD . About 40% of the 339 UDs we have de 
tected did not appear or disappear within the interval 
covered by the observations. Therefore, their lifetimes 
could not be computed. 

In the upper left panel of Figure |3] we show the position 
of appearance of all 339 UDs. The yellow circles indicate 
central UDs, the red squares peripheral UDs, and the 
blue triangles grain-origin UDs. Table [1] lists their mean 
properties. The average lifetime is about 18 min. This 
is relatively long co mpared with the values rep orted in 



previous works (e.g., Riethmiiller et al. 2QQ8b[ ), proba- 
bly because our manual procedure is capable of detect- 
ing fainter UDs. UDs move with an average speed of 
0.2-0.5 kms~^. The proper motion speed is defined as 
the distance between the points of appearance and dis- 
appearance divided by the time interval. Distances are 
not corrected for projection effects. The brightness ra- 
tio is the UD intensity (the maximum intensity within a 
±2 pixel area) relative to the intensity of the dark back- 
ground (/db)- We use the "dark background" — the region 
surrounding the UD but excluding the UD itself — as a lo- 
cal reference, /db is the mean intensity of pixels over a 
2" X 2" area centered in the UD whose intensity is darker 
than the average minus 0.5cr. Here, a represents the stan- 
dard deviation of the intensity within the 2" x 2" area. 
We find average brightness ratios from 1.5 (central UDs) 
to nearly 1.9 (grain-origin UDs). The mean UD diameter 
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Fig. 4. — Histograms of UD lifetime (upper left), average proper 
motion speed (upper right), brightness ratio (lower left), and di- 
ameter (lower right). Central UDs are displayed with black bars 
filled with yellow color, peripheral UDs with red color, and grain- 
origin UDs with blue color. The number of measurements in each 
category is shown in the upper-right corner of the histograms. 
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Fig. 5. — Intensity profiles emerging from the tip of a penum- 
bral grain (pluses) and from an UD (asterisks). The diamonds 
indicate the bisector positions at intensity levels of 10% and 23% 
(constituting x^high)? cind 49% and 62% (constituting ^'deep)• 

is 400-500 km. For each UD, the diameter is calculated 
as the average distance in eight radial directions along 
which the intensity is brighter than 1.2 /db^ oi" as the dis- 
tance to the cl osest inflection point. A similar method 
was adopted in [Watanabe et al. (2009a). There are some 
instances of zero diameter, which means that the UD in- 
tensity was darker than 1.2 /db- Brightness ratios and 
diameters are calculated for every time step within the 
UD lifetime, so in total we get 6279 individual values. 

Figure |4] shows histograms of these parameters for all 
the UDs detected in the observations. Central and pe- 
ripheral UDs are very similar except that the latter move 
faster. Grain-origin UDs, however, are clearly different: 
they have the fastest speeds, the highest intensity con- 
trasts, and the largest diameters. 

3.2. Derivation of Velocities: Line Bisectors 
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Fig. 6. — Original -Uhigh map observed at 08:31 UT (top) and 
subsonic filtered map (bottom). Negative velocities (blue) mean 
upflows along the line of sight, and positive velocities (red) mean 
downflows. The arrows mark the direction t o di sk center. The 
regions labeled A-C will be studied in Section 1X2] 



The mean wavelength position of the two line wings 
at a given intensity level is called the line bisector (Fig- 
ure [5| . Since different intensity levels sample different 
atmospheric layers, bisectors are often used to estimate 
the height variation of the line-of-sight velocity (e.g., 
[Tritschle r et al.||2QQ4[ [Ortiz et al^ |2Q10'). The 0% in- 
tensity level represents the line core, while 100% means 
the local continuum. 

We calculate bisectors only for Fe I 6301.5 A because 
the red wing of Fe I 6302.5 A is strongly blended with the 
telluric O2 6302.8 A line (see the Stokes / profile in Fig- 
ure [l]). The Fe I 6301.5 A line is less affected by blending, 
although its far blue wing (intensity levels > 70%) some- 



0.5 
0.4 

0.3 
0.2 
0.1 




O 

0) 
CD 

o 
CO 



0.06 
0.04 
0.02 
0.00 
-0.02 
-0.04 
-0.06 




-0.5 0.0 0.5 1.0 1.5 



-0.5 0.0 0.5 1.0 1.5 



0.04 

cn 

;^ 0.02 

^ 0.00 

CD 

B -0.02 

CO 

-0.04 




> 
cn 

CD 

o 
CO 



0.10 
0.05 
0.00 
-0.05 
-0.10 




-0.5 0.0 0.5 1.0 1.5 
wavelength - 6301.5 [A] 



-0.5 0.0 0.5 1.0 1.5 
wavelength - 6301 .5 [A] 



Fig. 7. — Stokes profiles observed at the position of an UD 
(marked with a triangle in Figure [ij and best-fit profiles returned 
by SIR (symbols and red solid lines, respectively). 



times show influence of molecular lines in the darkest 



umbral areas (Martinez Pillet & Vazquez 1993), which 
may result in a systematic blueshift. boi this reason 
we calculated four bisector positions at intensity levels 
of 10%, 23%, 49%, and 62%, avoiding bisectors close to 
the continuum. Each bisector position is obtained from 
a linear interpolation of the relevant intensities in the 
observed line profile. To reduce the noise, two bisector 
levels are averaged. As can be seen in Figure [5) this re- 
sults in two bisector velocities which will be called Vhigh 
(10% and 23%) and ^deep (49% and 62%). The Vhigh and 
v^QQp maps are not corrected for projection effects, but 
they mostly represent vertical flows because horizontal 
flows are usually weak in the solar photosphere. Nega- 
tive values of 'i^high and '^deep correspond to blueshifts (or 
upflows along the line of sight), and positive values to 
redshifts (or downflows). 

Following Title et al. ( 1989[ ), we apply a subsonic 
Fourier filter to the temporal sequence of 'Uhigh and 'i^deep 
to suppress disturbances with horizontal speeds larger 
than 4 kms~^, which are mostly due to p- modes and 
the residual temporal noise. In the umbra, a large 
scale ('^5000 km) velocity pattern corresponding to the 
p-mode oscillations is observed. This pattern is removed 
by the filter. We use an edge apodization of 10% in space 
and time to avoid the propagation of boundary errors. 
The effect of the subsonic filtering is demonstrated in 
Figure |6] The dominant redshift signal on the right side 
of the umbra at (x^y) ~ (37^^37'^) in the original im- 
age (upper panel) is absent in the filtered image (lower 
panel). However small local variations, which coincide 
with the UD's positions, keep their identity even after 
the subsonic filtering. 

For both I'deep and I'high, the zero velocity is deter- 
mined by averaging the filtered Vhigh map over umbral 
pixels with intensities below 0.4 /qg. This is done for 
each of the 68 frames of the sequence. The '^high and '^deep 
maps look similar, but the latter shows larger root-mean- 
square (rms) fluctuations. The typical rms velocities in 
the umbra are 0.11 kms~^ for '^high and 0.13 kms~^ for 
'Udeep- The rms variations increase significantly when the 
maps are not filtered. Thus, the filtering is essential to 
provide a good velocity reference, and also to remove the 
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large velocity offsets induced by the p-mode oscillations. 
This is particularly important when dealing with small 
velocities such as those observed in UDs. 

3.3. Derivation of Magnetic Properties: Stokes 
Inversion 

We determine the magnetic properties of the umbra 
by inverting the observed Stokes profiles with th e SIR 
code (Stokes Inversion based on Response functions; Ruiz 
Cobo fc del Toro Iniesta||1992D . The two lines are fitted 
simultaneously, excluding the telluric O2 blend in the red 
wing of the Fe I 6302.5 Stokes / profile. 

The inversion is carried out in terms of a one- 
component model atmosphere with constant (i.e., height- 
independent) magnetic fields and velocities, which is suf- 
ficient to explain the relatively symmetric Stokes profiles 
observed in the sunspot umbra (see Figure!?]). Zero stray- 
light contamination and unity magnetic filling factors are 
assumed. The inversion returns 9 free parameters: the 
three components of the vector magnetic field (strength, 
inclination, and azimuth), the line-of-sight velocity, the 
microturbulent velocity, and the temperature at 4 nodes. 
The initial guess model used t o start the inversion is the 
hot umbral model of Collados et al. I (0 ^994) . 

Sample maps of the retrieved magnetic parameters are 
displayed in Figure |3] The inclination and azimuth an- 
gles are expressed in the local reference frame to avoid 
projection effects. The inclination varies from 0° to 180° 
for vertical fields pointing away from and to the solar 
surface, respectively. The azimuth is measured coun- 
terclockwise from the positive x-axis of the figure. As 
expected, the sunspot shows an outward-directed radial 
magnetic field. We also note that the line-of-sight veloci- 
ties returned by the inversion agree well with the bisector 
velocities obtained from the 6301.5 A line. An animated 
version of Figure [3] is available in the on-line journal. 

4. CONVECTION IN THE UMBRA 

In the solar photosphere, the presence of strong mag- 
netic fields inhibits convective energy transport. This is 
why sunspots are dark and cool compared to the quiet 
Sun, which is covered by convective cells called granules. 
However, convection is not entirely suppressed even in 
the umbra: as a matter of fact, UDs are the manifesta- 
tion of a modified convective pattern. The properties of 
this pattern still need to be determined. Here we study 
how the characteristics of convection differ between the 
quiet Sun and the umbra, paying special attention to the 
morphology of the convective cells and the correlation 
between brightness and velocity. 

4.1. Morphology 

Figure [S] illustrates the morphological evolution of a 
central UD and a quiet-Sun granule. The images cover 
an area of 2^2 x 2^2. First we notice that the UD 
is much smaller than the granule (~0V3 vs IVO; cf., 
[Roudier fc Mullernl987( ). Second, UDs are relatively iso- 
lated whereas granules are closely packed between narrow 
(^0'.'3) inter granular lanes. Third, the intensity profile of 
granules is flat- top, sometimes with intensity depressions 
in the middle. A granule loses its identity (or disappears) 
by fragmentation or by merging with neighboring gran- 
ules. On the other hand, the intensity profile of UDs has 



a Gaussian shape, and they disappear mostly by fading 
out. 

Both granules and UDs exhibit a turbulent character: 
the structures displayed in Figure |8] change their bright- 
ness, barycenter, and shapes on timescales of only a few 
minutes. The variations are more pronounced in the case 
of granules. The average lifetime of UDs, 18 minutes, is 
slightly longer than the 5-15 minute duratio n of granules 
(e.g., [Bahng fc Schwarzschild 1961; Alissandrakis et al. 
1987[ lTitle et a 



I [Hirzberger et al.||19^| ~ 



4.2. Brightness vs Velocity 

The strong correlation between brightness and line- 
of-sight velocity in granular convection is well known. 
Bright areas (i.e., granules) show upflows, while dark ar- 
eas (the intergranular lanes outlining the granules) har- 
bor downflows. This is illustrated in the top panels of 
Figure [9] 

In the umbra, in addition to UDs, there exist diffuse 
areas with enhanced brightness. These areas may cor- 
respond to a convective pattern similar to (but weaker 
than) that of the quiet Sun, with UDs being another 
manifestation of the same pattern occurring at positions 
where convection is more vigorous. To test this possibil- 
ity, we have examined the correlation between brightness 
and velocity in the umbra: if bright structures are the 
result of convection, then they should preferentially be 
associated with upflows. The middle and bottom panels 
of Figure [9] show scatter plots for two umbral areas la- 
beled B and C in Figure [6j Area B is close to the light 
bridge and C represents the central umbra. In both areas, 
at 08:28 UT (left column) we find a tendency of upflows 
in bright regions and downflows in dark regions similar 
to that of the quiet Sun, yet with weaker correlation. 
For most of the scans with good seeing conditions, the 
slope k of the best linear fit to the velocity-brightness 
relation turns out to be negative, with average values 
of -0.11 ± 0.02 in area B and -0.24 ± 0.04 in area C. 
However, the negative correlation does not always per- 
sist in time and scans obtained a few minutes apart under 
good seeing conditions sometimes show opposite behav- 
iors (compare the lower panels of Figure [9]). For this 
reason, the present data do not allow us to unambigu- 
ously confirm the existence of a global convective pattern 
in the umbra. We will see later that the situation is dif- 
ferent for UDs. 

5. UMBRAL DOTS 

In this Section we give a detailed description of the 
evolution and statistical properties of individual UDs. 

5.1. Case Studies 

We select five UDs (UD#A-E) whose locations are in- 
dicated in Figure [lOj All of them were observed from 
appearance to disappearance. Movies of their temporal 
evolution ca n be f ound in the electronic journal. 

In Section [5. 1.5| another two UDs from one of the best 
scans of the sequence will be considered. We use them to 
demonstrate the existence of localized downflow patches 
around UDs in deep photospheric layers. 

5.1.1. Typical Central UD 

According to our visual inspection, more than 70% of 
the central UDs do not show flow field perturbations. 
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Temporal evolution of a central UD (top) and a quiet Sun gra nule (bottom). Displayed are continuum intensity maps for an 
area of 2'/ 2 x 2'/ 2. The coordinate system coincides with that of Figure ^ An intensity cut along the dashed line plotted at the center of 
the FOV is shown below each image. Times are given in UT in the intensity panels. 



i.e., no upflows or downflows are detected. Similarly, 
magnetic field perturbations associated with central UDs 
are usually not visible or very small. 17% of the central 
UDs split or merge with neighboring UDs. 
The temporal evolution of a typical central UD 



An example is the structure located next to UD#A at 
coordinates (34^1, 36^6), between T = 566 and 755 s. 
In Figure [l7, a local reduction of the field strength at 



(UD#A) is displayed in Figures 17 and 18 Figure 17 



shows maps of physical parameters starting 2 frames be 
fore the appearance of the UD (T=Os) and ending 2 
frames after its disappearance. In the left column of Fig- 
ure [18) curves representing the temporal variation of the 
parameters at the position of the UD (solid line) and the 
dark background (dashed line) are given. The curve for 
the dark background is an average over pixels with inten- 
sities lower than the average minus 0.5cr in the FOV of 
Figure The right column of Figure [18] shows spatial 
profiles along the x-direction at ^-positions co-moving 
with the UD. The Gaussian peak at x = 34!! 3 in the 
continuum plot corresponds to the UD. 

UD#A was born in a diffuse bright region, then in- 
creased in brightness (0s<T<377s), and finally merged 
with a neighboring UD (755 s<T<1006 s). As can be seen 
in both maps and plots, it did not show clear upfiows or 
downfiows. During most of its lifetime, the UD was lo- 
cated in a patch of enhanced redshifts whose morphology 
and amplitude do not correlate with the UD evolution. 
Other UDs in the FOV display localized upfiows of up 
to '^deep = 0-3 kms~^, but they do not persist in time. 



the position of UD#A can be seen during s<T<629 s. 
Then, the UD appears to gradually merge with a pre- 
existing weak field patch right above it. The amplitude 
of the magnetic field perturbations is very small, about 
50 G. The UD does not leave clear signatures in the in- 
clination maps. 

5.1.2. Distinct Central UD 



Figures 19 and [20| show the temporal evolution of a 
distinct central UD(UD#B). This UD differs from the 
others because of its large brightness and upflow. UD#B 
was born in a very dark area within the umbra and the 
continuum intensity became quite high about 9 minutes 
later (566 s<T<880 s). A second smaller intensity peak 
appeared in the interval 1069 s<T< 1321 s. Finally the 
UD faded in a diffuse bright background with no de- 
tectable intensity peak. There was a significant upfiow 
associated with the continuum intensity enhancement. 
The maximum upfiow of I'high ~ 0.3 kms~^ occurred at 
T=755s. A reduction of the field strength of the order 
of 50 G co-spatial with the UD can be observed during 
the first half of its lifetime. After T=755s, the region 
of weaker field strengths disappears and the UD collides 
with a pre-existing strong field region. Sometimes there 
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Fig. 9. — Bisector velocity -Uhigh vs continuum intensity for se- 
lected umbral and quiet Sun areas using the line scans taken at 
08:28 UT (left column) and 08:33 UT (right column). Negative ve- 
locities represent upflows along the line of sight. The positions of 
areas A-C are indicated in Figure [6] is the slope of the best 
linear fit and r the Pearson correlation coefficient. 

is a small spatial displacement (up to 0^2) between the 
brightness peak and the patch of reduced field strengths. 
A very small patch with more inclined fields appeared 
transiently from T = 566 to 755 s. 

The bright UD located at coordinates (31^5, 40^') also 
showed reduced field strengths and strong upflows (par- 
ticularly in ^high) for more than 6 minutes, from T = 
-125 to 251 s. 

5.1.3. Typical Peripheral UD 

Peripheral UDs are born in the peripheral region of the 
umbra, where the continuum intensity is brighter and the 
magnetic field more inclined. A significant property of 
peripheral UDs is their systematic motion toward the 
center of the umbra (see Table 1). Our visual inspection 
reveals that 55% of the detected peripheral UDs show 
upflows, but no systematic magnetic field perturbations. 
For example, 13% of peripheral UDs present a reduction 
of field strength, while 16% of them are associated with 
enhanced fields. Usually, these perturbations do not last 
more than a few minutes. 

Figures |2l] and [22] display the evolution of UD#C, a 
typical peripheral UD. The inward migration of this UD 
is very clearly seen in the online animation. The mi- 
gration speed is ~1.1 kms~^ during 0s<T<503s, and 
almost zero during 503 s<T<1635 s. In the final stages 
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Fig. 10. — Position of selected UDs (A-B: central UDs, C: pe- 
ripheral UD, D-E: grain-origin UDs). The background image shows 
the continuum intensity at 08:30 UT. The cross symbols indicate 
the trajectories of the UDs. The arrow marks the direction to disk 
center. 



of its life (1635 s<T<2390 s), the UD migrates again with 
a speed of 0.5 kms~^. The origin of UD#C is a diffuse 
bright area within the umbra. It shows three brightness 
peaks: at T = 377, 1132, and 1761 s. Only the first is as- 
sociated with strong upflows (I'high ^ —0.15 kms~^ and 
'^deep ^ —0.35 kms~^). The other two peaks show much 
smaller velocities, but still shifted to the blue compared 
with the dark background. A hint of redshifts can be 
seen in I'deep from T=2202 to 2390 s. 

The magnetic field shows a complex distribution. A 
patch of reduced field strength (~ —100 G) develops at 
the position of the UD and migrates inward to the umbra 
during 377 s<T< 1069 s. The magnetic field perturbation 
then disappears. At T=1635s the negative patch can be 
observed again at the position of the UD. It will survive 
for most of the remaining UD's evolution. In addition, 
another patch of increased field strengths is visible to 
the left of the UD, at {x,y) ^ (32'', 33''). This structure 
already exists from T=— 125 s (before the appearance of 
UD#C) and also migrates inward, but with slower speed 
(0.5 kms~^ during 0s<T<755s). The patch disappears 
when the UD's brightness decreases at around T=2202 s. 
The vector magnetic field of UD#C gets more inclined, 
especially during the migration. 

5.1.4. Typical Grain-Origin UDs 

Grain-origin UDs are characterized by inward migra- 
tion to the umbra center and high intensity contrasts 
(Table[T]). The perturbations associated with these UDs 
are more clearly visible than those of central and periph- 
eral UDs. Our analysis demonstrates that more than 
70% of the grain-origin UDs harbor upflows, while 40- 
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Fig. 11. — Close-up images of the central area of the umbra from the line scan taken at 08:30 UT. Clockwise, starting from top left: 
continuum intensity, -Uhigh? '^deep? field inclination, field strength, and bisector velocity at the 90% intensity level. Negative velocities 
represent upflows along the line of sight. The contours show continuum intensities of 0.4 /qs. 



50% show weaker and more inclined fields. The magnetic 
field perturbations occur preferentially near the penum- 
bra, becoming less prominent as the UD moves into the 
umbra. Sometimes one observes patches of both reduced 
and increased field strengths at the position of grain- 
origin UDs. 

In this section we describe the evolution of two typical 
grain-origin UDs: one from the disk center side of the 
spot and the other from the limb side. 

Figures [23] and [24] show the temporal evolution of the 
disk-center side UD (UD#D). It evolves from a filamen- 
tary structure into a circular shape and follows an un- 
usual trajectory: generally, grain-origin UDs move along 
the extension line of penumbral grains, but the trajectory 
of UD#D is almost perpendicular to it. The continuum 
intensity decreases a bit during 188s<T<314s, and then 
increases again. Upflows are observed in the first half 
of the UD evolution, weakening until they almost reach 
the background level at T = 692 s. The strongest flow 
('^deep< —0.6 kms~-^) occurs at the edge of the penum- 
bral grain early in the UD's evolution, around T = 62 s, 
and appears to be the continuation of the typical upflow 
of penumbral grains. 

UD#D is located at the boundary of weak and strong 
field regions. When it moves along the —y direction, the 
boundary also evolves as if the leading edge of the UD 
was always blocked by strong field walls (0s<T<817s). 
After that, the UD collides with the strong field region 
and then disappears together with the strong field walls. 
We could not detect any systematic perturbation of the 
field inclination caused by this UD. 

Figures 



the grain with an apparent speed of 0.7 kms~^. In the 
case of limb-side UDs, velocity perturbations are hard to 
detect. UD#E does show upflows in '^high and '^deep, but 
they are much weaker than those observed in center-side 
UDs. We attribute this to a line-of-sight effect working 
against the field-aligned flows in the inclined magnetic 
field. During the first part of its evolution, the UD is 
located in between two patches of weaker and stronger 
fields. Like in the case of UD#D, the migration seems 
to be impeded by strong field walls. Toward the end of 
the sequence (at T = 692 s), the UD only shows stronger 
fields than the surroundings. It does not seem to perturb 
the inclination of the umbral field. 

5.1.5. UDs with Downflow Patches 



UD#E, a 



25 and 26 show the temporal evolution of 



imb-side UD. This structure is detached from 



a penumbral grain and moves along the extension line of 



Ortiz et al. (2010) found evidence of downflow patches 
associated witn bright UDs in a pore. A bisector analysis 
of the Fe I 6301.5 line showed that these downflows are 
strongest in deep atmospheric layers. We tried to per- 
form a similar analysis, but this proved difficult for two 
reasons: 

1. In the darkest parts of our large spot, the Fe I 
6301.5 A bisectors at high intensity levels (> 70%) 
appear to show systematic blueshifts. 

2. Near the continuum, the quality of the bisector 
maps is very dependent on the seeing conditions. 

Therefore in this section we show examples of downflows 
extracted from the best scan in our data set, focussing on 
relatively bright areas around the center of the umbra. 

Figure [TT] shows a region of about 4" x 2" centered at 
(x, y) = (33^^ 38. 5'^). The two UDs marked in the upper 
left panel (UD#F and UD#G) exhibit no clear velocity 
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Fig. 12. — Histograms of average -Uhigh and -Udeep in central UDs 
(top), periplieral UDs (middle), and grain-origin UDs (bottom). 
The center-of-gravity value of each histogram is given in the upper 
right corner. 

signals in I'high and ^'deep^ but prominent upflows in the 
bisector map at the 90% intensity level (lower left panel). 
Next to those upflows, localized downflow patches can 
be seen in the vicinity of the UDs. The blueshift at 
the center of UD#F amounts to —1.3 kms~^, while the 
downflow patches have speeds of 0.16 kms~^ {-\-y side) 
and 0.75 kms~-^ {—y side). Similarly, the blueshift at the 
center of UD#G is —1.5 kms~^, and the downfiows to 
the top right attain 0.7 kms~^. The downflow patches 
have an approximate size of 0^2. In these UDs we do 
not see the central dark lanes predicted by [Schiissler fc 
iVoglerl (12006). 

UD^^F and UD#G show slightly weaker and more in- 
clined fields than their surroundings. When we inspect 
the continuum movie, these two UDs are both in the peak 
phase of their brightness. 

5.2. Statistical Properties 

In this section we describe the statistical properties of 
339 UDs (98 central, 112 peripheral, and 129 grain-origin 
UDs). Histograms of lifetimes, average proper motions, 
brightness ratios, and diameters have already been pre- 
sented in Figure [4j 

5.2.1. Bisector Velocities 

Figure [12] displays histograms of the mean I'high and 
Vdeep in UDs (averaged over an area of 7 x 7 pixels cen- 
tered at the UD's position), including all temporal steps 
within their lifetimes (i.e., 6279 samples). The center-of- 
gravity value is given in each panel. If UDs had no sys- 
tematic flows, the histograms would show a symmetric 
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Fig. 13. — Average value of -Uhigh (solid) and -Udeep (dash-dotted) 
as a function of radial distance from the UD center. Negative ve- 
locities mean upflows. The calculation is performed in the umbra, 
where the continuum intensity is darker than 0.4/qs. One pixel cor- 
responds to 0^06. The error bars show the la fluctuation within 1 
pixel bins. 

distribution about kms~^. This is the case for central 
UDs, with only a slight inclination to negative veloci- 
ties (upflows). On the other hand, a tail extending to 
strong upflows is seen in the histograms for peripheral 
and grain-origin UDs. The asymmetric distribution is 
more prominent in the case of grain-origin UDs. Sta- 
tistically, the upflows are stronger in fdeep compared to 
^high, indicating a deceleration with height. 

Figure [13] shows the average variation of the velocity as 
a function of radial distance from the center of the UD. 
Upflows decreasing outward are found in the region close 
to the UDs (<5 pixels), although the la fluctuation is 
large. The strongest upflows occur at the peak brightness 
position, with average velocities of 'yhigh= —0.07 kms~^ 
and Vdeep= —0.14 kms~^. Within a distance of 10 pixels 
(OV6) from the UD center, we always find upward veloc- 
ities on average, but no downfiows. 

5.2.2. Magnetic Parameters 

The local perturbations of field strength {AB) and field 
inclination (Ai) are obtained by subtracting a smoothed 
version of the maps from the original maps themselves. 
The smoothing is done with a boxcar of width 20 x 20 pix- 
els (870 X 870 km^), which is significantly larger than the 
typical UD size (see Figure |4|. 

Figure [14] shows the histograms of AB and Ai aver- 
aged over an area of 7 x 7 pixels centered at the UD's 
position. To our surprise, the histogram of AB indicates 
that the UD magnetic field is a bit stronger than the sur- 
ro undings, which seems to co ntradict the field-free model 
of Schiissler & Vogler (2006). The histo grams of Ai are 
almost symmetric, inaicatmg no preference for more ver- 
tical or more inclined fields at the position of the UDs. 

5.2.3. Temporal Evolution 

We have studied the temporal evolution of UDs us- 
ing the ones that were observed from birth to death and 
had lifetimes longer than 620 s. In total, 36 central, 50 
peripheral, and 66 grain-origin UDs were chosen for this 
analysis. By normalizing the lifetimes to unity it is possi- 
ble to average all the evolutionary curves. The results are 
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Fig. 14. — Histograms of AB and Ai for central (top), periplieral 
(middle), and grain-origin (bottom) UDs. These values are com- 
puted over an area of 7 x 7 pixels, centered at the position of the 
UD. Negative Ai means more vertical field lines, while positive Ai 
means more inclined fields. Shown in the upper right corner is the 
center of gravity of each histogram. 



shown in Figure 15 for eight parameters: peak brightness 
(/peak), brightness ratio, diameter, field strength, 'Uhigh, 
Vdeep^ and Ai. The parameters have been averaged 
over an area of 7 x 7 pixels centered at the UD's position. 

Evolution of Central and Peripheral UDs — For central and 
peripheral UDs, the peak brightness, the brightness ra- 
tio, and the diameter show a symmetric increase and 
decrease over time. The brighter /peak and weaker field 
strengths observed in peripheral UDs are a natural conse- 
quence of their location in the more external parts of the 
umbra. Both central and peripheral UDs harbor upfiows 
(negative 'i^high and '^deep) that grow with time, reach a 
maximum when the UDs are mature, and then decrease 
more or less symmetrically. The brightness follows a sim- 
ilar pattern. Thus, there is a positive correlation between 
brightness and upfiows, which is a sign of convection. 

The field strength and AB are nearly constant, al- 
though AB shows a tendency to increase during the evo- 
lution of peripheral UDs. Ai is around zero for central 
UDs and negative for peripheral UDs (indicating more 
vertical fields), with little variations over time. 

Evolution of Grain- Origin UDs — For grain-origin UDs, all 
parameters other than bisector velocities show patterns 
of monotonic increases or decreases. These patterns are 
caused by the smooth transition from penumbral grains 
to circular UDs in the central umbra. The bisector ve- 
locities peak shortly after the appearance of the UDs, 
when the tips of penumbral grains are completely de- 
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Fig. 15. — Temporal evolution of peak brightness (/peak) 5 bright- 
ness ratio, diameter, field strength, x^high? '^deep? AB^ and Ai for 
central (black), peripheral (red), and grain-origin (blue) UDs. The 
axis for grain-origin UDs is separately shown to the right in the 
panels of /peak 5 brightness ratio, and diameter. The standard er- 
rors of the mean (SEM) are indicated by the vertical lines. 



tached. Although with large scatter, locally weaker and 
more inclined fields are found in the first half of the UD's 
lifetime, while opposite properties (locally enhanced and 
more vertical fields) appear in the latter half. 

5.2.4. Scatter Relations 

Scatter plots of various parameters are displayed in 
Figure [T6| separately for central, peripheral, and grain- 
origin UDs (black, red, and blue, respectively). In gen- 
eral, all types of UDs exhibit similar relations between 
parameters, although there are some differences in bisec- 
tor velocities. Figure 16 shows that: 



1. Brighter UDs tend to have longer lifetimes. This 
is because most of them are recurrent. Other than 
that, the lifetime is almost constant regardless of 
the type of UD. 
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2. The diameter increases linearly for briglitness ra- 
tios up to 1.4 due to the way it is calculated. For 
larger ratios the diameter saturates to a constant 
value of about 500 km, demonstrating that UDs 
have a typical size. 

3. The proper motion speed decreases linearly with 
lifetime: short-lived UDs move faster. 

4. For peripheral and grain-origin UDs, the stronger 
upflows are found in the brighter structures. Cen- 
tral UDs show the opposite tendency. 

5. There is no correlation between AB and brightness 
ratio. 

6. Brighter UDs exhibit more inclined fields than 
dark UDs, although the inclination difference is not 
larger than 2° on average. 

7. There is a clear correlation between AB and field 
strength. For UDs with strong fields (>2000 G), 
AB is positive and reaches up to ^100 G. 

8. The most inclined fields are associated with neg- 
ative AB perturbations. UDs with more vertical 
fields tend to show positive AB values. 



Bisector velocities, AB^ and Ai are the most important 
parameters to decide on the convective nature of UDs. 
The field-free convection model of UDs suggests upfiows, 
weaker, and more inclined magnetic fields. These condi- 
tions are indeed true for relatively weak UDs (< 2000 G). 
However, strongly magnetized UDs (> 2000 G) usually 
show positive AB perturbations and more vertical field 
lines. We conjecture that the stronger upfiows observed 
in the darker central UDs may be the result of contami- 
nation by molecular lines. 

6. DISCUSSION AND CONCLUSION 

In this paper we have performed a detailed analy- 
sis of UDs in a mature sunspot using data from the 
CRISP spectropolarimeter. The excellent spatial reso- 
lution, temporal cadence, and polarimetric sensitivity of 
the measurements are ideal for UD studies. The pertur- 
bations caused by UDs are usually very small, and thus 
an statistical approach has to be followed to reveal their 
common properties. Our work addresses for the first time 
the temporal evolution of velocity and magnetic fields in 
and around UDs, using a statistically significant sample 
of UDs. 

Convection in the umbra — UDs are considered to be the 
manifestation of convection in the presence of the strong 
umbral field, while more vigorous convection occurs in 
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the quiet Sun in the form of granules. The morphologi- 
cal differences between granules and UDs can be seen in 
Figure [Sj Driven by overshooting cellular convection, 
granules are characte rized by sh arp edges and irregu- 
lar polygonal shapes fSpruit et al. 1990). UDs, on the 
other hand, show Gaussian brightness profiles. Linear 
theory reveals that the preferred horizontal scale of con- 
vection dec reases with increasing field strength (Weiss 
et al.||1990), which explains why UDs are smaller than 



granu. 



es. The convective origin of UDs seems well es- 
tablished, as many papers including ours found a good 
correlation between upflows and brightness. However, 
in the umbra there are other diffuse areas with enhanced 
brightness whose origin is still unknown. Our speculation 
that the diffuse bright areas of the umbra are also caused 
by convection could not be unambiguously confirmed on 
the basis of a unique relation between brightness and 
velocities (Figure [M. 



Photometric properties — IWatanabe et al. (2009a) found 
constant lifetimes regardless of the UD type and the 



structure of magnetic field at the position o f the U D. This 
is in agreement with our results (Section s. 2. 4| ). How- 
ever other studies report longer life times fo r brighter UDs 
( [Tritschler fc Schmidt, 2002^ Bharti et al pOlOJ , and we 
speculate this is partly due to the fact that bright UDs 
tend to reoccur at the same position. 

A correlation between shorter lifetime and faster 
proper motion is reported for the first time in this paper 
(Figure 16). If the energy dissipation rate is proportional 
to the UD speed, the lifetime can be expected to be re- 
duced for fast-moving UDs, as observed. We also find 
that the travelled distance depends linearly on lifetime, 
i.e., long-lived UDs travel longer distances even though 
they move at lower speeds. The same conclusion can 
be obtained from a sim ilar analysis of the data set of 
Watanabe et al | ( |2009aD . 

We found UD diameters consistent with the values re- 
ported in the literature, i.e., about 400km on average. 
The Gaussian shape of the histogram (Figure |4| and the 
lack of dependence of the diameter on brightness ratio 
(Figure 16) suggest that UDs indeed have a "typical" 
size, regardless of their type. This common UD size is 
probably determined by a universal near-surface strati- 
fication in mature sunspots. Howev er, the scatter plot 
analysis performed in Section 5.2.4 did not reveal any 
physical parameter having a strong correlation with the 
UD size. 

The in t ensity oscillations in UDs reported by, e.g.. 



Rimmele 



(1997) have not been studied in this paper 
because of the uncertainties that residual seeing fluc- 
tuations may introduce. However it is tru e that many 
UDs s how recurrence, as observed also by [Louis et al.| 
([2012"). For example, the peripheral UD#C displayed 
in Figure 22 reappeared twice within a time interval of 
13 minutes. This timescale is co mparable to the oscil- 
latory period of the UD s shown in Rimmele (1997) and 
Watanabe et all (|2009b|). 



Categorization of UDs — We classified the observed UDs 
in central, peripheral, and grain-origin UDs according to 
their place of birth. Do these categories represent phys- 
ically different structures or different manifestations of 
the same phenomenon? Grain-origin UDs have larger 
brightness ratios, larger sizes, and faster proper motions 



than the other UDs. The temporal evolution of grain- 
origin UDs is smooth and shows monotonic changes, 
while central and peripheral UDs show mound-shaped 
evolutionary curves (Figure 15). Despite this, the scatter 
plots presented in Figure [16] suggest that the properties of 
grain-origin UDs lie on the extension lines of those of cen- 
tral and peripheral UDs. The differences between them 
may arise from stronger convection in grain-origi n UDs 
rendered p ossible by the weaker background field. |Spruit| 
& Scharmer (2006) speculated that field- free convection 
can explain both U Ds and penumbra l grain s. T he com 



puter simula tions of Heinemann et al. (2007) and Rempel 
et al. (2009*) succeeded in reproducing basic properties o: 



penumbral filaments and UDs as weakly-magnetized con- 
vective structures. Our results are in general agreement 
with the predictions of this scenario, although they also 
indicate that UDs are far from being completely field- 
free (at least in the photospheric layers accessible to the 
observations) . 

Substructures — Localized downflow patches at the pe- 
riphery of UDs are considere d to be a signature of over - 
turning convection in UDs (Schiissler & Vogler 2006). 
In Section [5'.1.5| we presented some examples of down- 
flow patches without performing a full statistical anal- 
ysis. The patches of Figure 11 have sizes of 0^2 and 
redshifts of up to 0.75 kms^V^Both the size and the 
velocity are in go od agreement with those reported by 
Ortiz et al. (2010) in a pore. The fact that the downflow 
patches are observable only at very high bisector levels 
(i.e., deep photospheric layers) is also consistent with the 
results of those authors. 

However, many UDs do not show downflow patches in 
our data. This lack of detection could be due to: 

1. Insufficient spatial resolution. 

2. The existence of downflows only in deep layers that 
cannot be probed by the Fe I 6301 and 6302 A lines. 

3. The transient nature of the downflows, which could 
appear only in a particular phase of the UD's evo- 
lution. 

4. The possibility that the convective energy escapes 
to the upper layers instead of returning to deep 
layers (see the narrow jet-like u pflows above the 



cusp in Schiissler fc VQgler||2006 ) . 



The two UDs featured in Figure 11 are in the peak 
brightness phase when they show downflows. Possibly 
the speed of the downflows reaches a maximum when 
the brightness is also maximum. 

Velocities in UDs — Upflows and brightness follow similar 
evolutionary patterns in UDs (Figure fl5]). T his correla- 
tion support s the convective natur e ofUDs (Sobotka fc 
Jurcak 2009) [ Watanabe et al .|2010|). The scatter plots of 
velocity vs brightness shown in Figure [16] also point to a 
convective origin of UDs. For central UDs, however, we 
observe stronger blueshifts in darker structures, although 
the tendency is not very pronounced. We suspect this is 
an artifact caused by systematic blueshifts in very cold 
umbral areas. 

The velocities observed within UDs are likely to repre- 
sent field- aligned flows, because upflows are readily found 
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on the disk-center side where the field is closer to the 
line -of- sig ht direction, compared to the limb side (Sec- 



tion 



5.1.4). The same effect can be observed in the Ev- 



ershed How (Figure [3|, which is also a field-aligned fiow. 
Magnetic field in and around UDs — The field-fre e convec- 



tion model of UDs (Schiissler & Vogler 2006) predicts 
weaker and more inclined fields in UDs. Our scatter 



analysis (Figure 16) confirms these properties, but only 
for UDs with fieldsDelow 2000 G. In strongly magnetized 
UDs (> 2000 G), the magnetic field is enhanced and more 
vertical compared to the surroundings. For grain-origin 
UDs, the physical conditions also depend on the phase 
of evolution: in the first half of their lifetime, weaker 
and more inclined fields appear, while stronger and more 
vertical fields are observed in the latter half as the UDs 
intrude into the umbra. To the best of our knowledge, 
the enhanced and more vertical fields of strongly mag- 
netized UDs cannot be explained by currently available 
UD models. 

We observe strong field regions at the migration front 
of grain-origin UDs for the first time (see the evolution of 
UD#D and UD#E in Section [5T4| . These strong field 
regions seem to impede the migration of the UDs. The 
situation is remin iscent of that modeled by [SchlicheiT 
maier et al. (1998), where a weakly magnetized penum- 
bral tlux tube pushes and compresses the pre-existing 
vertical field at the leading edge. However, also a weakly 
magnetized convective structure would produce a com- 
pression of the adjacent magnetic field, which has to wrap 
around th e field- free gas. The M HD simulations per- 
formed bylHeinemann et al. ( 2007[ ) predict the existence 
of enhanced field regions surrounding grain-origin UD 
only in layers deeper than the continuum forming region 
(see Figure 3 in that paper), but on both the leading 
and the tail sides. Our observation did not find field 
enhancements on the tail side. 



Final remarks — This work extends our knowledge of the 
temporal evolution of velocities and magnetic fields in 
UDs. We found some new and unanswered results that 
may provide constraints to future modeling efforts. A 
pioneering comparison of observational and computer- 
simula ted UDs has been performed by Bharti et al.] 
(2010), and this kind of studies should be extended. 

At the same time, more spectropolarimetric observa- 
tions of UDs at high cadence should be performed. The 
temporal resolution of our data, 63 s, seems appropri- 
ate to track the evolution of UDs, but is insufficient for 
resolving the evolution of UD substructures. UDs will 
remain one of the most challenging targets for solar ob- 
servations in the coming years. 
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Fig. 17. — Temporal evolution of a central UD (#A) from 2 frames before appearance until 2 frames after disappearance. From top to 
bottom, maps of continuum intensity, ^'high5 '^deep? field strength, and field inclination are shown. T = s corresponds to the appearance 
of the UD. The black plus symbols indicate the UD position. The {x,y) coordinates follow the coordinate system used in Figure [l] 
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Fig. 19. Same as Figure [T7| for central UD#B. 
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Fig. 20. — Same as Figure [Ts] for central UD# B. The times corresponding to the five different fines in the right panels are 503 s (black 
plus, dashed line), 629s (black square), 755s (red diamond), 880s (blue square), and 1069s (blue plus, dashed line). 
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Fig. 21. — Same as Figure for peripheral UD#C 
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Fig. 22. — Same as Figure [TS] for peripheral UD#C. The times corresponding to the five different hnes in the right panels are 817s 
(black plus, dashed line), 943s (black square), 1132s (red diamond), 1384s (blue square), and 1510s (blue plus, dashed line). 
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Fig. 23. 



Same as Figure [TtI for grain-origin UD#D 
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Fig. 24. — Same as Figure [TS] for grain-origin UD#D. The times corresponding to the five different fines in the right panels are 62 s 
(black plus, dashed line), 251 s (black square), 503s (red diamond), 692s (blue square), and 880s (blue plus, dashed line). 
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Fig. 25.— 



Same as Figure 1171 for grain-origin UD^^E. 
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Fig. 26. — Same as Figure [TS] for grain-origin UD#E. The times corresponding to the five different fines in the right panels are 62 s 
(black plus, dashed line), 251 s (black square), 503s (red diamond), 692s (blue square), and 880s (blue plus, dashed line). 



